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f^^ ' We discuss finite temperature quantum Monte Carlo methods in the framework of 

^\1 ' the interacting nuclear shell model. The methods are based on a representation of 

the imaginary-time many-body propagator as a superposition of one-body prop- 
agators describing non-interacting fermions moving in fluctuating auxiliary fields. 
yJ , Fermionic Monte Carlo calculations have been limited by a "sign" problem. A 

CaJ ■ practical solution in the nuclear case enables realistic calculations in much larger 

configuration spaces than can be solved by conventional methods. Good-sign in- 
teractions can be constructed for realistic estimates of certain nuclear properties. 
We present various applications of the methods for calculating collective properties 
and level densities. 
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^D • 1 Introduction 

o 
o\ 

C^ I A variety of models have been developed over the years to explain the observed 

^^ . properties of nuclei in various mass regions. One of the more fundamental of these 

^P ' models is the interacting shell model, where valence nucleons (outside closed shells) 

(^ . move in a mean field potential and interact via a residual nuclear force. This 

effective interaction can be traced back to the bare nucleon-nucleon interaction via 
Brueckner's G-matrix. i^ i^ 

The nuclear shell modeLwas applied successfully in the description of jB, saa, 
513 . and lower p/-shell nucleiBoLl. However, the size of the model space increases com- 

binatorially with the number of valence nucleons and/or orbitals, and conventional 
diagonalization of the nuclear Hamiltonian in a full shell is limited to nuclei up to 

At finite temperature many levels contribute, and diagonalization methods are 
even more limited. However, there are methods that do not require direct diago- 
nalization of the Hamiltonian. At temperature T the nucleus is described by its 
free energy F{T) = -TlnZ{T), where Z(T) = Twr"/^ is the partition func- 
tion expressed in terms of the nuclear Hamiltonian H . It is difficult to calculate 
the partition function because H includes interactions that are strong, and non- 
perturbative methods are required. Thermal mean-field approximations, e.g., the 
Hartree-Fock approximation, are useful and tractable. Small amplitude fluctua- 
tions around the mean field can also be taken into account in the random phase 
approximation (RPA). However, in the finite nuclear system, large amplitude fiuc- 
tuations are important. In fact, interaction effects can be taken into account exactly 
when all fluctuations of the mean fleld are properly included. This is formally ex- 
pressed by the Hubbard-Stratonovich transformation^ In this transformation the 
nuclear propagator in imaginary time is represented as a superposition of one-body 
propagators that describe non-interacting fermions moving in fluctuating auxiliary 
fields. Quantum Monte Carlo methods were developed to carry out the integration 






X 



Tenth International Conference on Many-Body Theories 



over the large number of auxiliary fields .Dfl These calculations are computation- 
ally intensive and became feasible only with the introduction of parallel conmuters. 
Similar methods were developed for strongly correlated electron systems.clll3 

The initial applications of the Monte Carlo techniques were severely limited 
by the "sign" problem, which is generic to all fermionic Monte Carlo methods 
and occurs for all realistic niiclear interactions. We have developed a practical 
solution to this sign problemB that made possible realistic calculations in much 
larger configuration spaces than could be treated previously. 

Other quantum Monte Carlo methods were developed for the many-nucleon sys- 
tem. Variational and Green function Monte Carlo methods were used [Successfully 
to derive properties of light nuclei from the bare nuclcon-nucleon forceJlj Recently, 
a hybrid Monte Carlo method was suggested, in which the spin variables are de- 
coupled using the Hubbard-Stratonovich transformation^ 

The shell model Monte Carlo (SMMC) methods are reviewed in Section g. Some 
of their applications are presented in Sections ^, ^ and ||. 

2 Methods 

2.1 The Hubbard-Stratonovich transformation 

The shell model Hamiltonian contains a one-body part and a residual two-body 
interaction. A generic Hamiltonian with two-body interactions can be brought to 
the form 



H = '^eana + -^Vapl , (1) 



2 

a a 

where ea is the single-particle energy in orbital a, pa are linear combinations of 
one-body densities pij — ajoj, and Va are interaction "eigenvalues." 

The canonical density operator e~^^ , where f3 — 1/kT is the inverse tempera- 
ture, can also be viewed as the many-body evolution operator in imaginary time /3. 
The Hubbard-Stratonovich (HS) transformatioro is an exact representation of this 
propagator as a path integral over one-body propagators in fiuctuating external 
fields (j{t). It is derived by dividing the time interval (0, /3) into Nt time slices of 
length A/3 each, and rewriting e~^^ = (e""^^^) '. The propagator for each time 
slice can be written as e"-^^^ « Ha e^'^'^"""" Ha e-^^^^"'^' to order (A/5)2. Each 
factor with p^ can be written as an integral over an auxiliary variable CTq 



where Sq. = ±1 for Va < and Sa — H for Va > 0. By introducing a different set 
of fields aa{Tn) at each time slice r„ = nAf3, we obtain the HS representation of 
the propagator 



e-/3ff^ 



Jv[a]G^U^, (3) 



where 



G. = e-^^^S=."'''°I'^°(-^ (4) 



Tenth International Conference on Many-Body Theories 



is a Gaussian weight, and 

U„ = e-^'3'''(^"' . . . e-^'3'''(^i' = Texp | - / h„{T)dT I (5) 



/3 



is the propagator of non-interacting nucleons moving in time-dependent external 
one-body auxihary fields CTq (t) ■ The one-body Hamiltonian in (0) is given by 



hair) = ^eafla +^Sa\Va\aa{T)pa , 
a a 

and the metric in the functional integral (0) is 



I?[CT]=n„,„ da^{Tn)V^l3K\/2TT 



(6) 



(7) 



The thermal expectation of an observable O can be represented from the HS 
transformation (pi) as 



(O) 



Tr(0e-'3«) _ JV[a]G„{0)^TrU^ 



Trfe-/3ff) 



JV[a]G„TrU„ 



(8) 



where {0)„ = Tr (0[/cr)/Tr C/o- is the expectation value of O for non-interacting 
particles moving in external fields crir). 

The calculation of the integrands in (ph can be reduced to matrix algebra in the 
single-particle space. To show that, we denote by TVs the number of valence single- 
particle orbits, and represent the one-body propagator C/^ in the single-particle 
space by an N^ x Ns matrix \J^. The grand-canonical trace of Ua in the many- 
particle fermionic space can be expressed directly in terms of \J^ 



TrUa =det(H-U<, 



(9) 



Eq. (|9|) is just the grand-canonical partition function of non-interacting fermions 
(moving in time-dependent external fields a). 

Similarly the grand-canonical expectation value of a one-body operator O = 
Si Oija\aj can be calculated from 



{a\aj)c 



U„ 



(10) 



The grand-canonical expectation value of a two-body operator can be calculated 
using Wick's theorem {alar^aiat;) = {ajai) (a'^ak) ~ {a\ak){ajai) together with (|lO|). 
In practice, we are interested in canonical expectation values, for which the 
number of particles A is fixed. Canonical quantities can be calculated by an exact 
particle-number projection. The canonical (one-body) partition function is given 
as a Fourier transfornJlj 



Tr^f/, = 



J-/3M 

N, 



E 



-i<timA 



det [l ■ 



P«0m „/5m 



u.] , 



(11) 
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where (jj^n = I'Km/Ns {in = 0, . . . , Ng) are quadrature points and /i is a chemical 
potential. Similarly for a one-body observable O 

—f^fiA ^s /I 

Tr^ iOU.) = ^ E --""-^ det [l + e^^-^^^V^] tr ( , ^ ^..,,.,,^-1 O^ 

(12) 
where O is the matrix with elements Oij. In practical calculations it is necessary 
to project on both neutron number N and proton number Z , and in the following 
yl wiU denote {N,Z). 

2.2 Monte Carlo methods 

The integrands in (H) are easily calculated by matrix algebra in the single-particle 
space (see, e.g., Eq. (ph and (p^)). However, the number of integration variables 
(Ta{T) is very large. For small but finite A/?, this multi-dimensional integral can be 
evaluated exactly (up to a statistical error) by Monte Carlo methods.Ll 
The first step is to define a positive-definite weight function 

W, = G„\TtU^\. (13) 

We can rewrite the expectation value (||) of an observable O as 

where 

$, = TrC/,/|Tr{/,| (15) 

is the sign of the one-body partition function. 

In the Monte Carlo approach, a random walk is performed in the space of 
auxiliary fields a = {(Ta{Tm)} that samples the cr-fields according to the distribution 
Wa-. The M^- weighted average of any cr-dependent quantity Xc^ can be estimated 
from 

^^^^"^ ^ jD[a]W, - M 1.^-""' ' (16) 

where ct*^'^-' are M uncorrelated samples. The statistical error of {Xcr)w can be 
estimated from the variance of the "measurements" X^ik) . Though a standard 
random walk can be constructed by the Metropolis algorithm, a modification based 
on Gaussian quadratures improves the efficiency in the nuclear case. til 
Using (^ and (p^), thermal expectation values are calculated from 

(0) = M^ . £|i%i. . (17) 

Often the sign $0- in dl3) is not positive for some samples a, and we must 
"measure" it together with the observable (O)^- 
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2.3 The Monte Carlo sign problem 

When the average sign {^a)w is smaller than its uncertainty, the Monte Carlo 
method fails. This leads to the so-called Monte Carlo "sign" problem: the sign of 
the integrand fluctuates among samples, and the integral is the result of a delicate 
cancellation that cannot be reproduced with a finite number of samples. Often the 
average sign decreases exponentially with /3, and the problem becomes more severe 
at low temperatures. The sign problem is generic to all fermionic Monte Carlo 
methodsu and has severely limited their applications. In particular, the problem 
occurs for all realistic effective interactions in the nuclear shell model. We have 
developed a practical solution to the sign problem in the nuclear case that allows 
realistic calculations in very large model spaces.Q 

Certain interactions can be shown to have a "good" sign. A time-reversal in- 
variant Hamiltonian with two-body interactions can be written as 

a a 

where pa is the time-reverse of pa and Va are real. When all Va are negative in the 



representation (18), the grand-canonical partition function is positive-definite for 
any sample a and the interaction has a good sign. The proof is based on a general- 
ization of Kramer's degeneracy to the complex plane. The one-body Hamiltonian 
that appears in the HS decomposition for ( p^ ) has the form 



^^^<^afla+^{VaSaPa+VaSaPa) ■ (19) 



When all Wq < 0, Sq = 1, and the Hamiltonian (|l^) is time-reversal invariant 
ha- = h„. The eigenstates of Uq- appear then in time-reversed pairs with complex 
conjugate eigenvalues {Ai, A*}, and the grand-canonical partition function Tr U^ = 
det(l -I- Ucr) = Hi|l -I- Aip is positive-definite. The canonical partition function for 
even-even nuclei is also positive-definite. 

Realistic effective nuclear interactions have both positive and negative "eigen- 
values" Va, leading to a severe sign problem at low T. However, the eigenvalues 
that are large in magnitude are all negative and therefore have a good sign. This 
property of realistic nuclear forces can be used in several ways to overcome the sign 
problem. 

In general, the Hamiltonian is decomposed into "good" and "bad" parts: H — 
Hg + Hb- The good Hamiltonian contains the one-body part plus the two-body 
terms with Va < 0, while the bad Hamiltonian contains the two-body terms with 
Va > Q- We construct a family of Hamiltonians 

Hg=HG+ gHB (20) 

that depends on a continuous coupling parameter g. For any g < 0, the Hamilto- 
nian (|o|) is good-sign and accurate Monte Carlo calculations can be performed for 
{0)g = Tt {Oe-'^"^)/Tr (e-f^""). The dependence of {0)g on g is expected to be 
"smooth" and can therefore be extrapolated to g = 1. We found, through tests in 
the sd- and lower p/-shell, that low-degree polynomial extrapolations (usually first 
or second order) are sufficient. The extrapolation technique is used in Section pi 
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Certain nuclear properties such as collectivity and level densities can be reliably 
calculated-hy constructing nuclear Hamiltonians that are completely free of the sign 
problem .113113 These Hamiltonians includ£_correctly the dominating collective com- 
ponents of realistic effective interactions.EZI An example of a good-sign Hamiltonian, 
used in the applications discussed in Sections || and ||, is 

H^Y. '-f''^ - 9°P^°''^^ ■ ^^°''^ - ^ E ^^0(^'°) • 0(^'°) , (21) 



where 



(A:T) 



ah 

o''^^' - -^j^^ll^jA'^y^UM. X a,j(^^-) . (22) 



^j.rn.mt 



The modified annihilation operator is defined by 
(^_p-ni+2-'"ta^ _^ _^^^ and a similar definition is used for p(-^'^). The 
single-particle energies. Eq are calculated from a Woods-Saxon potential plus 
spin-orbit interaction.E3 V in (Uw is the central part of this potential, and the 
multipole interaction is obtained by expanding the surface-peaked interaction 
w(r, r') = --x{dV/dr){dV/dr')6{r — f'). Since the monopole pairing and isoscalar 
multipole- multipole interactions are all attractive, they lead to a good-sign 
Hamiltonian (|2l|). 

Experimental gap parameters (found from odd-even mass differences) are used 
to determine the pairing strength go through a particle-projected BCS calculation. 
The parameter x of th^ surface-peak interaction is determined self-consistently: 
X^^ — Jr. dr r^ {dV/dr) [dp/dr), where p is the nuclear density. The renormaliza- 
tion factors k\ account for core polarization. 

2.4 Particle-number reprojection 

The SMMC methods are computationally intensive. Calculations can be done more 
efficiently in the particle-reprojection techniqueJlH Suppose we use Monte Carlo 
methods to sample a nucleus with A particles. The ratio Z^i/Zj[ between the 
partition function of a nucleus with A' particles and the partition function of the 
original nucleus can be calculated from 



Za{P) Tr^e-/3« \ Tr aU^ 

where (. . .)w is defined in (pq). Each of the partition functions within the brackets 
in ( p3| ) is calculated from Q. The expectation value of an observable O for the 
nucleus with A' particles is calculated using 

Tt.,OU„\ lTr.,U, 



)(^ 



Ir 4/ U (T 
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In the particle-reprojection method, the Monte Carlo saniphng is carried out by 
projecting on a nucleus with fixed A, and the partition functions and observables 
are calculated for a family of nuclei with A' j^ A using Eqs. (|2|) and (p4|). 

In the following Sections we present some of the applications of the shell model 
Monte Carlo approach. 



3 Quenching of the Gamow- Teller strength 

Gamow- Teller rates, especially in the iron xegion, are important inputs in models 
of stellar collapse and supernova formation.E3 They determine the electron-capture 
rates, which in turn affect the electron-to-baryon ratio. 
The Gamow- Teller transition operator is given by 

GT^ oc CTr± , (25) 

where a and r are the spin and isospin operators, respectively. The GT~^ transition 
operator changes a proton into a neutron, and its strength function can be measured 
in {n,p) reactions.L^ The total strength B{GT^) is found experimentally to be 
strongly quenched relative to estimates based on the independent-particle model. 
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Figure 1. Total Gamow- Teller strength B{GT^) in p/-shcll nuclei. The solid syrphols are the 
renormalized SMMC strengths and the open symbols are the experimental results. EH After Ref. 
23. 



Calculations thaj—include up to 2p — 2/i correlations account for some of the 
observed quenching.E3 To estimate how much of the quenching can be attributed 
to proton-neutron correlations in the ground state, we perfomied SMMC calcula- 
tions in the complete p/-shell for nuclei in the iron region.ES The total strength 
is renormalized by (1/1.26)^ ~ 63% to account for the normalization of the axial 
coupling constant from its free value oi qa — 1-26 io qa — ^ in the nuclear medium. 
The renormalized strengths are shown in Fig. ^ in comparison with the experi- 
mental values. We see good agreement across the shell. We remark that a similar 
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renormalization of the Gamow- Teller operator in complete sd-shell calculations also 
leads to an agreement with the data.El 

We conclude that any observed quenching beyond an overall renormalization of 
~ 63% is due to proton-neutron correlations. This has been confirmed in conven- 
tional shell model calculations in the p/-shell, where the Gamow- Teller strength is 
found to decrease as a functimi of the number of particle-hole excitations that are 
included in the model spaceO 

In SMMC it is also possible to calculate imaginary-time response functions and 
invert them by the maxent method to obtain strength functions.l3 Such methods 
were used to calculate Gamow- Teller strength distributions for p/-shell nucleiHa 

4 7-soft nuclei 
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Figure 2. Free energy surfaces of ^^''Xe and ^^*Te at several temperatures. The surfaces are in- 
ferred from shape distributions calculated in SMMC. The inset on the right shows typical intrinsic 
shapes sampled in the Monte Carlo together with their standard deviation. After Ref. 15. 



One of the outstanding problems in nuclear structure is to demonstrate that 
various types of collectivity (known from phenomenological models of nuclear struc- 
ture) emerge in microscopic calculations. In particular, nuclei in the 100 - 140 mass 
region are known to be 7-soft, i.e., their potential energy is insensitive with respect 
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to non- axial deformations. This region corresponds to the 50 - 82 shell (for both 
protons and neutrons), and could not be solved previously in fully microscopic cal- 
culations because of the prohibitively large size of the model space. Using SMMC 
we have provided the first microscopic evidence of softness in this mass regionll3. 
We have used the good-sign Hamiltonian ( pi| ) and, among the various multipolc in- 
teractions, kept only an attractive quadrupole interaction renormalized to account 
for core polarization. A quadrupole pairing interaction (in addition to monopole 
pairing) was necessary to obtain correctly the excitation energy of the first 2^ state 
in tin isotopes. 

To determine the quadrupole shape distribution P((7^), we calculated for each 
sample a the average (Q^)^ and variance A^ = {Q'^)a — {Q)a of the quadrupole 
mass operator Q^ = J2^^y2ti- It is convenient to convert the shape distribution in 
the intrinsic frame_P((7g, q'2) to an effective free energy surface at finite temperature 
using the relation^ P{q'„,q'2)dq'Qdq'2 ex exp[-F{P,r,T)/T](3'^\sm{S-/)\dpdj, where 
/3 and 7 are the intrinsic quadrupole shape parameters. Such free energy surfaces 
are shown in Fig. ^ for ^^^Te and ^^^Xe at several temperatures. The inset shows 
typical intrinsic shapes, demonstrating the weak dependence of the free energy on 
7- 

5 Level densities 

Level densities are important for theoretical estimates of nuclear reaction rates. In 
particular, level densities are necessary input for nucleosynthesis calculations. The 
s and r processes are determined by the competition between neutron capture and 
beta decayp] and the neutron-capturaxates are proportional to the level density of 
the corresponding compound nucleus .E3 i_. 

Experimental data on level densities are available from different methods:Ej di- 
rect counting of levels at low energies, neutron and proton resonance data and 
charged particle spectra at intermediate energies, and Ericson's fluctuations at 
higher energies. 

Most conventional calculations of level densities are based on Fermi gas models 
in which important correlations are neglected. It is often found however, that 
empirical modifications of the Fermi gas formula (also known as Bethe's formulaEj) 
can describe well experimental data. Particularly useful is the backshifted Bethe 
formula (BBF) where the ground state energy is backshifted by an amount A. The 
backshift parameter originates in pairing correlaticpa and shell effects. A modified 
version of the BBF is due to Lang and Le CouteuiEil 



/TT 



p{E,) ^ 9^a~HE, - A + tyie^V-iE.-^) ^ (26) 

where < is a thermodynamic temperature defined by E^ — A = at^ — t and g — 2. 
The experimental level densities for many nuclei are well described by (EQ) when 
both a and A are adjusted for each nucleus. It is therefore difficult to predict the 
level density to an accuracy better than an order of magnitude. 

The nuclear shell model provides a suitable framework for the microscopic cal- 
culation of level densities. The dimensionality of the model space is often too large 
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to allow for exact diagonalization, and truncations that might be appropriate for 
low-lying states are not suitable at finite excitation energies. Instead we use the 
Monte Carlo methods described in Section g. 

The calculations of level densities in SMMCE3 are based on a thermodynamic 
approach. The level density is obtained from the canonical partition function by 
an inverse Laplace transform 

piE) = r ^eP'^ZiP) . (27) 

J -too 27ri 

In practice we are interested in the average level density. It can be obtained from 
( p7[ ) in the saddle point approximation 

p{E) « (2^/3-2c)-i/2eS ^ (28) 

where the canonical entropy S and heat capacity C are given by 

S = \nZ + l3E; C = 0\{H^) - {Hf) = -fi^dE/dfi . (29) 

In SMMC we calculate the thermal energy E{0) as an observable {E{(3) = 
{H)fj), and then find the partition function by integrating the exact thermodynamic 
relation E = —d\nZ/d(3: 

ln[Z(/3)/Z(0)] = - / dl3'E{(3'). (30) 







Z(0) = Tr 1 is the total number of states in the model space. The entropy and 
heat capacity are then calculated from Eqs. (|2^ ) to give the level density (|2^). 

5.1 Level densities in the iron region 

We have used the methods of Section |^ to calculate level densities for nuclei in the 
iron region. The model space includes the pf + gg /2-shell, and is good for excitation 
energies up to E^ ~ 20 MeV. With the inclusion of the 5(9/2 orbit we can calculate 
both even and odd parity states. 

We used the good-sign Hamiltonian (^ and kept only the quadrupole, octupole 
and hexadecupole terms with renormalization factors of ^2 = 2, ^3 = 1.5 and 
^4 = 1. Using experimental odd-even mass differences we determined go ~ 0.2 
MeV. 

Fig. H demonstrates the calculation of the level density of ^^Fe. The top and 
bottom insets show the entropy versus energy and the heat capacity versus /3, re- 
spectively. For comparison we show the Hartree-Fock (HF) results by dashed lines. 
In the HF we observe a discontinuity in the heat capacity around /3 ~ 1.3 MeV~^ 
due to a shape transition from a deformed to spherical nucleus. This discontinuity 
is washed out by the fluctuations in the SMMC calculations. The total SMMC level 
density for ^^Fe is shown in Fig. versus excitation energy Ex (solid squares). The 
solid line is inferred from the experiment .E3 We see good agreement between the 
microscopic calculations and the data, without any adjustable parameters. 

The SMMC level density can be well-fitted to the BBF. We can therefore extract 
a and A from the microscopic calculations. Fig. shows the parameters a (left) 
and A (right) as a function of mass number for even-even nuclei in the mass range 
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Figure 3. Total level density of ^®Fe versus excitation energy E:^'- SMMC (solid squares), exper- 
iment (solid line) and HF (dashed line). The top inset shows the entropy S{E), and the bottom 
inset is the heat capacity C{0). The solid squares are the SMMC results while the dashed lines 
are obtained in the HF approximation. After Ref. 16. 



A '^ 50 — 70.e3 The SMMC results (solid squares) are compared with the empirical 
values of Ref. 34 (solid lines). We observe that a varies smoothly with mass, but 
that A shows shell effects and is enhanced at A^ = 28 or Z = 28. 
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Figure 4. Single-particle level density parameter a (left) and backshift parameter A (right) as a 
function of mass number A, obtained by fitting the SMMC level densities to the BBF. The solid 
lines are the empirical values of Ref. 34. Shown are results for even-even nuclei from Fe to Ge. 
After Ref. 33. 
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5.2 Level densities by particle-number reprojection 

The projection on an odd number of particles introduces a new sign problem, even 
for good-sign Hamiltonians. Moreover, the calculations of level densities are time- 
consuming since each nucleus requires new Monte Carlo calculations for all tem- 
peratures. 



The particle-number reprojection method of Section 2.4 can be used to carry 



out the Monte Carlo sampling for an even-even (or N ^ Z) nucleus A followed by a 
reprojection to find the thermal energies for a series of nuclei A'. We have applied 
this method to calculate the level densities of manganese, iron and cobalt isotopes 
by reprojecting from ^^Fe and ^*Co.tj Since the Hamiltonian (pi| ) depends on A 
(mostly through its single-particle spectrum), suitable corrections should be made 
to the thermal energy. 
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Figure 5. Level densities of cobalt isotopes versus excitation energies obtaijied in the particle- 
number reprojection method. The solid lines are fits to the BBF formula (|26|). After Ref. 19. 



Results for the level densities of cobalt isotopes, reprojected from ^^Fe, are 
shown in Fig. |^. We observe that the level density of an odd-odd cobalt (e.g., 
^■^Co) is higher than the level density of the subsequent even-odd cobalt (^^Co) 
although the latter has a larger mass. The systematics of a and A is shown in Fig. 
|6| for manganese, iron and cobalt isotopes, including odd-A and odd-odd nuclei. 
The SMMC results (solid squares) are compared with the experimental values (x's) 
and the empirical formulas of Ref. 34 (solid lines). The staggering effect in A 
originates in pairing correlations. The microscopic calculations are often in better 
agreement with the data than are the empirical values. 
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Figure 6. The single-particle level density parameter a (left) and backshift parameter A (right) 
in manganese, iron and cobalt isotopes. The values determined by paa-Wcle-number reprojection 
SMMC (solid squares) are compared with the experimental values (x's) E2I and the empirical values 
of Ref. 34 (solid lines). After Ref. 19. 



5.3 Parity-projected level densities 

We have calculated the parity dependence of level densities by introducing par- 
ity projection in the HS representation. The projection operators on even- and 
odd-parity states are given by P± = (1 ± P)/2, where P is the parity operator. 
The parity-projected thermal energies E±{f3) = Tr(i?P±e"'^^)/Tr(F±e"^^) can 
be written as 



E±{f3) 






(31) 



where C(a) = Tr U^, Cp(ct) = Tr (PUa) and {H)p^ = Tr {HPUa)/Tr {PU„). PU, 
can be represented in the single-particle space by the matrix PUo-, where P is a 
diagonal matrix with elements (— )^'' (^a is the orbital angular momentum of the 
single-particle orbit a). Consequently, the integrand in ( pi| ) can be calculated by 
matrix algebra in the single-particle space. From E±{f3) we calculate the positive- 
and negative-parity level densities using the method discussed in Section 0. 

The left panel of Fig. R shows the parity-projected level densities for ^^Fe. 
We see that p+ 7^ p_ even at the neutron resonance energies. This is contrary 
to the common assumption of equal even and odd parity densities, often used in 
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nucleosynthesis calculations. 

It is difficult to calculate the odd-parity level density for the even-even nuclei at 
low excitation energies because of a new sign problem introduced by the projection 
P_. The sign problem affects less the odd-to-even partition function ratio. This 
ratio is calculated from 






1- 






w. 



1 



Cp(a) 



and shown in the right panel of Fig. [7| for several nuclei in the iron region. 



(32) 




10 15 20 



E^ [MeV] 



|3(MeV'') 



Figure 7. Left: parity projected level densities p± in ^^Fe. Right: odd-to-even partition function 
ratio versus /3 for ^®Fe, ^'^Ni arid ^^Zn. The symbols are the SMMC calculations and the solid 
lines are obtained from Eq. (B3) of the statistical quasi-particle model (see text). After Rcfs. 16 
and 35. 



The parity distribution as a function of temperature (or excitation energy) can 
be understood in terms of a simple statistical model.E3 The single-particle states are 
divided into two groups of odd- and even-parity states. Assuming that the particles 
occupy the single-particle states randomly and independently, the probability to 
find n nucleons in states with parity tt (vr is the parity of the group with the 
smaller occupation) is a Poisson distribution, P(n) = /"e~-^/n!, where / = (n). 
For an even-yl nucleus, an even (odd) parity many-particle state is obtained when 
n is even (odd) . We find a simple formula for the odd-to-even parity ratio 



tanh/ 



(33) 



Z+m E„oven^(") 

/ is estimated from the Fermi-Dirac occupations of deformed single-particle states 
with parity tt. The distribution P(n) can be calculated in SMMC by projection 
on states with parity tiH. In the iron region we find good agreement with the 
Poisson distribution above the pairing transition temperature of '^ 1 MeV. Below 
this-temperature, our model is still valid but for the quasi-particles. / is now given 
byN / = Eae. /a = Eae.Il + e'^''"]"' where K = ^{t^ - A)^ + A^ are the quasi- 
particle energies (and A is the chemical potential) . The solid lines in Fig. M are the 
result (|3^) of the quasi-particle model. 
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Conclusions 

We have discussed a quantum Monte Carlo approach to the interacting shell model 
at finite and zero temperature. The initial applications of these techniques have 
been severely limited by the Monte Carlo sign problem. A practical solution to 
the sign problem in the nuclear case has allowed new realistic calculations that 
could not be carried out in conventional diagonalization methods. Alternatively, 
good-sign interactions can be constructed for realistic estimates of certain nuclear 
properties (e.g., collective properties and level densities). We have presented a 
variety of applications: the quenching of the Gamow- Teller strength, microscopic 
evidence of 7-softness, and accurate calculations of level densities and their parity 
distribution. 

Finally, it is important to keep in mind that the Monte Carlo approach is com- 
plementary to the conventional shell model approach; it cannot be used to calculate 
detailed spectra but on the other hand it is very useful for finite and zero temper- 
ature calculations. 
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